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Abstract 

We argue that two-step models, like String and Glasma models, with creation first 
of sources extended in rapidity that after locally decay into particles, lead to long 
range forward-backward correlations due to fluctuations in the number or the colour 
of the sources, and to a similar expression for the F-B correlation parameter b. 

In the simplest String percolation scenario, b increases with centrality or the 
number of participating nucleons. However, asymptotically it decreases with energy. 
This is different from the Glasma model where b increases with both, centrality and 
energy. 
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The study of particle correlations and fluctuations can give an important con- 
tribution to the understanding of the behaviour of matter at high density in 
high energy heavy ion collisions. Interesting experimental results have been 
recentely presented, [1,2,3,4], as well as various theoretical developments (see 
[5] and references therein). In this paper we discuss forward-backward rapid- 
ity correlations in the general framework of models with particle production 
occuring in a two-step scenario. 

In models where particle production occurs in two steps, formation first of a 
rapidity structure reflecting the longitudinal distribution of the colour fields 
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(as it happens with the Glasma [6] or with the dual string model [7]), followed 
by local emission of particles, the presence of long distance rapidity forward- 
backward correlations is, in a sense, unavoidable. Even if particle emission 
from the created sources is totally uncorrelated, as usually assumed, long range 
correlations occur in the process of averaging over fluctuations in the colour or 
number of the emitting sources. In multi-collision processes, as nucleus-nucleus 
collisions, such fluctuations are naturally expected [8,9]. 

In the simplest model, all the extended in rapidity sources are equal, and 
fluctuate in their number or colour N, and n particles are locally produced, 
with the result [9]: 

(n) = Nn , (1) 

where (n) is the average multiplicity, and n the single source average multi- 
plicity, and 

r,2n ^2 N^-N 2 Id 2 

B2 /<«) 2 = -^ + ^-, (2) 

where D 2 = (n 2 ) — (n) 2 and d 2 = n 2 — n 2 . The first term in the right hand side 
of (2) represents the rapidity long range correlation, and the second term the 
local short-range one. It is clear that the dynamics associated to the initial 
conditions of the collision is controled by the first (long range correlation) 
term in (2). The second term controls local production (resonances, short 
range correlations, etc.). For simplicity we shall assume that the local, single 
source particle distribution is a Poisson, which means 

d 2 /n = 1 . (3) 



From (2) and (3) we obtain 



D 2 = ^ + (n), (4) 



K 
with 

TV 2 

K = _ _ 2 . (5) 
N 2 -N 

It should be noticed that the quantities n, (n),N and K are, in principle, 
dependent on the width and position of the rapidity window under study. We 
shall concentrate on symmetrical, AA, collisions. 

If we consider a window in the forward direction (or backward direction) we 
obtain 

(n f ) = (n b ) = Nn , (6) 

and 

D)f = ^ + <»/> • (7) 
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If we consider next two symmetrically placed windows, one forward and one 
backward, separeted by a rapidity gap, such that the short range contribution 
is practically eliminated, we have 



Dl, = ^f- (8) 



The correlation parameter b, in the relation (n^p = a + bnp, is given by 

b = D 2 bf /D 2 ff , (9) 

and we obtain 

t = ITW' (10) 

This result was previously obtained in two different models, [10] and [11], 
having in common the presence of sources extending along the rapidity axis, 
and fluctuations in their number or overall colour. 

In the Glasma approach, [11], the long range correlation piece behaves as l/a s , 
while the short range correlation one behaves as a s , such that K/(nf) ~ 
a 2 , a s being the QCD coupling constant. The F-B correlation parameter b 
increases with centrality or the number of participating nucleons, N part , the 
same happenning with the increase of the energy. 

In the String approach we do not have such straightforward result: one has to 
analize the quantities (rif) and K. This is what we will do next. 

It is well known that dual string models cannot reproduce correctely the height 
of the rapitity plateau [12] and the N part dependence of dn/dy [13] without 
inclusion of string fusion, percolation, or any other multiplicity suppression 
mechanism. 

In string percolation, the relevant parameter is the transverse density of strings, 
V, 



^=(l) 2jv (v^), (11) 

where r is the string radius, R the transverse radius of the overlap region and 
N(y/s) the number of strings at a given energy. Making use of the standard 
relations in high energy interactions with nuclei, [12,13], 



N A 4/3 /N A 1/3 



N~N p (^) ,R^R p ^) ^ (12) 
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Fig. 1. D^j, Eq.(7), and Dfj, Eq.(8), as functions of N part , by making use of (16) 
and (18). Data are from [3]. Errors were estimated from the figures of [3]. 

where N p is the number of strings in pp collisions and R p the nucleon radius, 
we obtain 



(13) 



For the moment, as we are only interested in the N part dependence, we write, 
from (13), 

V = AN 2 p ' a l t , (14) 

A being a parameter. From recent and not so recent work [14], it is known 
that for Npart — 350, central Au-Au collisions, rj ~ 3, which means A ~ 0.06. 
We fix this value for A. 

The most important consequence of percolation is the presence of the colour 
reduction factor F(rj), [15], 



F( V ) 



'1 - e - " 



(15) 



such that (1), for the forward region, can be rewritten in the form 

(n f ) = F(AN 2 p £ t )BN*l 3 rt , (16) 



B being a free parameter. Note that asymptotically F(rj) — > r] 1 ^ 2 . 
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Fig. 2. a) (up) as a function of the number of participants N part . The curve is from 
(16) B = 0.119 ±0.006. 

Regarding the inverse of the normalized variance of the source distribution, 
K, it is known to increase with N part , for N part >30 [2,4]. We simply assume 
that string correlations, for large density, are proportional to the density i] : 
N ~ 7i, (W - N 2 - N) ~ 77, such that 

K~r], (17) 

and we write, see (14), 

K = CN 2 ^ t , (18) 



C being a new parameter. 

Note that relations (17) and (18) in the framework of percolation should not 
apply in the small 77 region: the K{rj) should have a minimum and even increase 
with i], as i] — > 0, as seen in pp collisions [16]. We can take that behaviour into 
account, but it does not change the final result as RHIC data are not sensitive 
to small N part . 

Making use of (18) and (16) with (15) one sees that, for large N part , 

K/(n f )^N£i\ (19) 

which means that the correlation parameter b increases with centrality, or 
the number of participants, as in [11]. In the no-percolation limit, F(r)) = 1, 
K/ ( n f) ~^ NpaJt ', which means that b shows a faster increase with centrality. 
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Fig. 3. K as a function of N part . The curve is from (18) with C = 10.1 ± 1.3. 

Recentely, the STAR collaboration at RHIC, [3], has presented results for Djj- 
and Dip shown in Fig.l, as a function of N part in Au — Au at y/s = 200GeV. 
Making use of the experimental points and of equations (7), for D 2 j f) and (8) 
for D%f, we have extracted the "experimental" points for (n/) (Fig. 2) and K 
(Fig. 3), and fitted them with equations (16), for (rif), and (18), for K (see 
the curves in Figs. 1,2 and 3). Note that K(N part ) in Fig. 3 is not very much 
constrained by data, in particular in the low N part region, as mentioned above. 

It is interesting to ask another question: what happens to the F-B correlations 
at fixed N part , with increasing energy? The answer is, see (13), 



(n f ) = F(A'N p (^))B'N p (V^) , 



(20) 



and 



K = C'NpiyG) , 



(21) 



A', B' and C being parameters, such that, asymptotically, 

K/(n f ) — > iV p (v^) 1/2 , (22) 

which means that the F-B correlation parameter b decreases as y/s goes to 
infinity. In the absence of percolation, F{rj) = 1, b is essentially constant. 
In the Glasma model, [11], b increases with the energy. We may thus have 
here a possibility of distinguish the Glasma model from the String Percolation 
model. As in moving from RHIC to LHC the number of strings in supposed 
to increase by a factor around 2 (see, for instance, the second reference in 
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[14]), we expect at LHC a reduction of b, for central collisions and similar 
phase-space, of the order of 25% in comparison with RHIC results. However it 
should be mentioned that this result strongly depends on the validity of (17) 
(see [16]). If K ~ r] a , with a < 1/2, then b also increases with energy in the 
percolation scenario. 
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